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ABSTRACT 

We have discovered a concentration of extremely red objects (EROs; R — K > &) m 
the field of the z = 2.69 quasar QSO 1213-0017 (UM 485), which is significantly over- 
abundant compared to the field ERO surface density. The optical/near-IR colors of the 
EROs and numerous other red galaxies in this field are consistent with elliptical galaxies 
at z = 1 — 2. HST optical images for a subset of galaxies show regular morphologies, 
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most of them being disky or diffuse and without any obvious evidence for interactions. 

Ground-based IR images show similar morphologies, indicating any dust reddening in 
these objects is spatially uniform. Optical spectroscopy with the W. M. Keck Telescope 
has found that four of the red galaxies lie at z ^ 1.31, and a fifth lies in the foreground 
at 2; = 1.20. Of the z ^ 1.31 galaxies, one is a reddened AGN while the remaining 
three have rest-frame UV absorption-line spectra characteristic of old (few Gyr) stellar 
populations, similar to the old red galaxy LBDS 53W091 at z = 1.55. Including the 
Mg II absorber seen in the QSO spectrum, we find five galaxies at z ~ 1.31 spread over 
1.5 h'^Q Mpc on the sky. These results suggest we have discovered a coherent structure 
of old galaxies at high-redshift, possibly associated with a massive galaxy cluster. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: clusters: general 
— galaxies: evolution — galaxies: stellar content — infrared: galaxies 

1. Introduction 

The observational study of galaxy formation and evolution can broadly be divided into two 
approaches: (1) searches for primordial galaxies, typically oriented toward high rcdshifts, in or- 
der to scrutinize formation and evolution processes as they occur, and (2) studies of existing old 
galaxies to decipher their origin and life history from their current properties. With the advent 
of infrared imaging detectors, a population of infrared-bright, extremely red objects (EROs) have 
been uncovered which is relevant to both of these approaches. 

Elston et al. (1988) discovered a few objects with exceptional optical/IR colors {R — K>5)m 
the first deep near-IR sky survey. While the brightest objects were shown to be z < 0.8 ellipticals 
(Lilly et al. 1988; Elston et al. 1989), their tantalizing conjecture that the fainter {K > 18) red 
objects are z > 1 ellipticals has remained unresolved until recently. Subsequent deep infrared 
imaging of primarily high-redshift radio galaxy and quasar fields serendipitously uncovered handfuls 
of objects with extreme colors (McCarthy et al. 1992; Eisenhardt & Dickinson 1992; Graham et al. 
1994; Hu & Ridgway 1994; Soifer et al. 1994; Dey et al. 1995; Djorgovski et al. 1995). Deep 
multicolor optical/IR sky surveys have begun to assemble larger samples of these objects using 
wcU-dcfincd selection criteria, an essential starting point for statistical studies (Cowic et al. 1994; 
Hall & Green 1998; Cohen et al. 1999; Barger et al. 1999; McCracken et al. 2000). In particular, 
the recent availability of 1024 x 1024 pixel IR detectors has made possible surveys for significant 
numbers of EROs (Thompson et al. 1999). 

The use of the designation "extremely red" has varied in the literature. This was especially 
true in earlier work where the description was applied to any object appearing in IR images but 
undetected optically. In this paper, we follow Graham &; Dey (1996) and use the term "extremely 
red object" (ERO) for a source with R — K>Q. This criteria was an operational one, encompassing 
several IR-detected, optically-invisible sources known at the time, and it has since become widely 
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used. In particular, Thompson et al. (1999), who have conducted the widest deep R — K survey to 
date, also adopt this definition. They confirm EROs defined by this criteria are unusual objects, 
being the reddest 2% of the K < field galaxy population. However, we will also pay attention 
in this paper to galaxies with R — K > 5, the criteria used by Elston et al. (1988) and Cohen et al. 
(1999), as this is about the expected color of passively-evolving elliptical galaxies at z > 1. We will 
generically refer to this larger sample as "red galaxies." 

EROs are very optically faint (R > 24) which has hampered studies of these objects. Successful 
spectroscopy to determine their redshifts and physical nature has only become possible with the 
development of the Keck 10-m telescope. Recent work has found some members of this "missing 
population" are luminous (> L*) galaxies at 2; > 1; hence determining their origin has direct 
relevance to the formation of massive galaxies and AGN. Being solely defined by optical/IR color, 
the few EROs with measured redshifts form a heterogeneous population as expected, comprising at 
least two broad classes: (1) ultraluminous dusty star-forming systems, perhaps akin to local objects 
like Arp 220, and (2) massive galaxies with old passively-evolving stellar populations. 

The best studied ERO to date, ERO J164502-K4626.4 (object 10 of Hu &: Ridgway 1994, 
hereinafter "HRIO"), is the prototype for star-forming EROs. This object is a 2; = 1.44 dusty, 
ultraluminous galaxy with enormous ongoing star formation (> 1000 Mq yr~^) suggested by its 
sub-mm continuum emission (Graham & Dey 1996; Cimatti et al. 1998; Dey et al. 1999). However, 
since HR 10 has among the reddest colors (I — K = 6) of known EROs, it is unrepresentative of the 
bulk of the population or at least is an extreme example. These EROs may provide excellent case 
studies of optically obscured star formation in the early Universe, especially given the association 
of faint EROs with some of the sub-mm emitting sources found by SCUBA (Small et al. 1999). 
In fact, it may be that intense, very brief bursts of star formation are a common mode of star 
formation at high redshift. 

There are also examples of EROs as galaxies with old stellar populations. Passively evolving 
ellipticals at z > 1 arc expected to have large R — K colors, which stem from their rising spectral 
energy distributions (SEDs) longward of 4000 A being redshifted into the near-IR. Therefore, 
selection using very red colors is an excellent method to search for early- type galaxies at 2; > 1. 
This is especially true in searching for clusters since the surface density of galaxies becomes quite 
high at faint magnitudes and it would otherwise be hard to distinguish a cluster from the foreground 
and background populations. However, ellipticals at these distances are expected to be optically 
very faint so measuring spectroscopic redshifts is challenging. Dickinson (1995) has identified a 
cluster of elliptical galaxies associated with the powerful radio galaxy 3C 324 at z = 1.206 along 
with a foreground structure at z = 1.15. To date, the highest-redshift collection of old, red galaxies 
spectroscopically confirmed has been found by Stanford et al. (1997). They have discovered a 
z = 1.27 cluster by its large J — K colors; the cluster galaxies have R — K >^ and rest-frame UV 
spectra resembling local elliptical galaxies. There is also a neighboring cluster at z = 1.26 found 
from its X-ray emission by Rosati et al. (1999) which contains spectroscopically old red galaxies. 
Isolated examples of old EROs have been discovered at still higher redshifts. The very weak radio 



-4- 



sources LBDS 53W091 at z = 1.552 (Dunlop et al. 1996; Spinrad et al. 1997) and LDBS 53W069 
at z = 1.432 (Dunlop 1998; Dey et al. 2000) both have R - K ^ 6 and rest-frame UV spectra 
which imply ages of a few Gyr. Recently, Soifer et al. (1999) have identified an R — K 7 object 
as an old galaxy at z = 1.58 based on associating a large continuum break observed at Ril jim with 
redshifted 4000 A break. 

Discovery oi z > 1 galaxies with old stellar populations offers several powerful lines of in- 
quiry into understanding galaxy formation and evolution and its cosmological context. Detailed 
comparison of the absorption lines and continuum breaks of these galaxies with galaxies at z = 
may prove fruitful in tracing the evolutionary course and enrichment history of the oldest stellar 
populations. Absolute age dating of these galaxies would provide a constraint on the time scale of 
galaxy formation and the age of the Universe. Moreover, clusters of old galaxies at high rcdshift 
can provide testing grounds for competing scenarios of galaxy formation; the predicted appearance 
of early- type galaxies in these clusters is dramatically different in hierarchical galaxy formation 
scenarios as compared to monolithic collapse ones (e.g., Kauffmann &; Chariot 1998). Finally, the 
existence of these old galaxies at high redshift potentially can constrain cosmological parameters 
and theories of structure formation (e.g., Peacock et al. 1998). 

We are conducting an on-going study of the nature of these EROs, using deep optical and 
near-IR imaging from the ground to assemble a large sample of EROs for statistical study. We 
have been acquiring high-resolution morphological information from HST optical and Keck near-IR 
imaging, and we are obtaining spectroscopy from Keck in the optical and near-IR to determine ERO 
redshifts and physical properties. Deciphering the identity of EROs based on comparing broad- 
band colors alone to theoretical stellar population synthesis models is dubious given that the model 
parameter space is vast, at least comprising age, metallicity, and reddening variations. Spectroscopy 
is essential. In addition, given the apparent heterogeneity of the population, a reasonably large 
sample of objects needs to be studied to understand the nature and relative abundances of the 
subsets, instead of the spectroscopy of individual EROs which has been done to date. 

In this paper, we present a study of the EROs in the field of QSO 1213-0017 (RA = 12^15™49.8^ 
Dec = -00°34'34" ; J2000.0). This z = 2.69 quasar, also known as UM 485, has exceptionally strong 
and complex Mg II absorption systems at z = 1.3196 and z = 1.5534 (Steidel Sz Sargent 1992). In 
§ 2, we describe optical and near-IR imaging covering an 11 arcmin"^ region around this field and 
Keck optical spectroscopic follow-up of galaxies selected by their very red colors. We consider in 
§ 3 the surface density of the red galaxies, their morphologies, and their spectroscopic redshifts. 
We examine in § 4 the collection of red galaxies as a whole, both their spatial distribution to con- 
sider the possibility that they are members of a cluster at z = 1.31 and their spectrophotometric 
properties to understand their stellar populations. We summarize our findings in § 5 and discuss 
their implications. 

Throughout this paper, we assume a cosmology with = 1, A = 0, and Hq = 50 h^o km s~^ Mpc 
At z = 1.31 with these parameters, 1" = 8.58 h^Q kpc; the luminosity distance (II = 9.40 h'^Q Gpc; 
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and the angular diameter distance = 1.78 h^Q Gpc. 

2. Observations 

2.1. Optical Imaging 

Optical imaging of the field of QSO 1213-0017 was obtained in April 1993 using the KPNO 
4-m Mayall telescope as part of a program to image z > 1 Mg II absorbing galaxies near QSO 
sightlines (see Steidel 1995). We used the PFCCD equipped with a 2048 x 2048 Tektronix CCD 
and a pixel scale of O'.'47/pixel. A total integration of 3500 s was obtained through the 7^-band 
filter (Ac = 6930 A, AA = 1500 A) under photometric conditions and l'.'25 FWHM seeing. The 
data were reduced using standard techniques and calibrated onto the AB magnitude system^ using 
spectrophotometric standard stars from the lists of Massey et al. (1988) and Oke & Gunn (1983). 
The TZ filter is a compromise between the standard Cousins Rc and Ic filters. A filter trace for 
TZ and a photometric transformation to the Kron-Cousins system are given in Steidel & Hamilton 
(1993). For the very red objects considered in this paper, TZ ~ Rc after accounting for the color 
terms, so the TZ — K colors given for the EROs in this paper are roughly equivalent to Rc — K 
colors given for other EROs in the literature (see § 3.1). 

We obtained four images of the 1213-0017 field with the WFPC2 camera aboard Hubble 
Space Telescope (HST) on U.T. 1998 January 11. Data were taken using the FSIAW filter with 
a total integration of 4500 s. The integrations were taken in slightly offset positions, allowing 
reasonable identification of cosmic ray hits in the individual images which were then excluded from 
the averaging to create the final image. 

2.2. Near-IR Imaging 
2.2.1. KPNO 

Infrared data were obtained in February 1994 on the KPNO 4m Mayall telescope using the 
IRIM imager, equipped with a 256 x 256 NICMOS-3 HgCdTe detector and sampling the sky at 
C.'GO/pixel. A total integration time of 3660 s was obtained in a non-repetitive grid, in which 
the telescope was moved after each 60 s of integration time. Conditions were photometric with 
I'.'O FWHM seeing. The data were reduced using the DIMSUM package (Eisenhardt, Stanford, 



^Note that our preliminary results in Liu et al. (1999) designated the TZ filter as "Rs" and presented Vega-based 
magnitudes, instead of AB ones. For the 72, filter, the offset from AB mags to Vega-based mags is -0.28 mag, i.e.. 
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Dickinson, & Ford, priv. communication) in IRAF.^ The data were calibrated using the UKIRT 
faint standards (CasaU & Hawardcn 1992), resulting in Vega-based magnitudes. The data were 
taken in the Ks filter (McLeod et al. 1995) which is slightly bluer and narrower than the standard 
K filter. As discussed by Hall et al. (1998), the shift between Kg and Kukirt is expected to be 
quite small, less than 0.04 mag for objects as red as H — K = 1. In the absence of IR colors for all 
the objects to compute color terms, we assume Ks = K hereinafter, which will have a negligible 
effect on any of our results. 

2.2.2. Keck 

We obtained deeper, higher spatial resolution IR imaging of two sub-fields on U.T. 1998 
May 14-15 using the facility near-IR camera NIRC (Matthews &: Soifer 1994) of the 10-meter 
W. M. Keck I Telescope located on Mauna Kea, Hawaii. The camera employs a Santa Barbara 
Research Corporation 256 x 256 InSb array and has a pixel scale of O'.'150/pixel resulting in a 
38" field. We observed 2 fields close to the quasar, one to the southwest (containing R4, R5, and 
R6) and one to the northeast (containing R7, R9, and RIO). The SW field, centered about (-27" 
E, -7" N) from QSO 1213-0017, was observed with the standard JHK filters with integration 
times of 900, 2040, and 2280 s, respectively, totaled over the two nights. The NE field, centered 
about (27" E, 10" N) from QSO 1213-0017, was observed in J for 780 s and in K for 1200 s on 
15 May 1998 UT. We used integration times of 6 or 20 s per coadd, depending on the filter. After 
one minute of integration was coadded, the sum was saved as an image, and then the telescope was 
offset by a few arcseconds. The telescope was stepped through a non-redundant dither pattern, 
and an off-axis CCD camera was used to guide the telescope during the observations. Conditions 
were non-photometric while observing the NE field and also for part of the SW field observations. 
For the SW field, we culled frames which were non-photometric from the reduction process. For 
the NE field, the data are used only to examine morphologies. 

The data were reduced in a manner typical for near-IR images. An average dark frame was 
subtracted from each image to remove the bias level. We constructed a fiat field by median averaging 
scaled images of the twilight sky. A preliminary sky subtraction of the images was performed to 
identify astronomical objects. Then for each image, we subtracted a local sky frame constructed 
from the average of prior and subsequent images, excluding any of the identified astronomical 
objects from the averaging. We used the brightest source in each field to register the reduced 
frames, which were then shifted by integer pixel offsets and averaged to assemble a mosaic of the 
field. Bad pixels were identified by intercomparing the registered individual images and masked 
during the construction of the mosaic. All these reductions were done using custom software written 
for Research System Incorporated's IDL software package (version 5.1). 



IRAF is distributed by the National Optical Astronomy Observatories, which arc operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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We observed the faint HST IR standard SJ 9143 (Persson et al. 1998) as a flux calibrator so 

the resulting magnitudes are Vega-based. Aperture photometry of the registration object in each 
individual frame verified the SW field data which were retained were nearly photometric, with any 
systematic errors of order 5% or less. There are no obvious point sources in the images, but the 
most compact objects have a FWHM of (X'5. 

2.3. Optical/Near-IR Photometry 

To identify very red objects, we compiled a photometric catalog of all objects in the KPNO 
K-hand mosaic using the SExtractor software of Bertin & Arnouts (1996), version 2.0.21. Each 
pixel in the iiT-band image was multiplied by the square root of its exposure time to create a 
mosaic with uniform noise over the entire field. Objects in the noise- normalized mosaic were then 
identified as any set of contiguous pixels 1. 5a above the background level with an area equal 
to a FWHM-diameter circular aperture (10 pixels), i.e., a highly significant detection, and then 
aperture photometry was done on the original K-hand image. We used the resulting "MAG_BEST" 
magnitudes, which for uncrowded objects use apertures determined from the moments of each 
object's light distribution. This method is similar to that of Kron (1980) and is designed to recover 
most of an object's flux while keeping errors low. The TZ image was transformed to be registered 
with the K image, and photometry was done using the same apertures for each object as for the 
iC-band image. The Galactic reddening towards this field is small, E{B — V) = 0.02 (Schlegel et al. 
1998), so we apply no extinction corrections to the photometry. 

We ran Monte Carlo simulations to verify the accuracy of the measured magnitudes and errors. 
We inserted 10 artificial stars of known magnitude and color into the images and then processed the 
images with SExtractor in the same fashion as the original data. For each magnitude and color bin, 
we ran the simulation 30 times for a total of 300 stars and then compared the recovered magnitudes 
and colors to the input ones in order to determine the random and systematic photometric errors. 
The artificial stars were chosen to span a range in magnitude and color encompassing all the real 
objects, with 0.5 mag steps in K magnitude and 1 mag steps in TZ—K color. A low-order polynomial 
surface was then fitted to the photometric errors from the simulations as a function of K magnitude 
and TZ — K color, and the fit coefficients were used to assign photometric errors for the real objects. 
In all cases the systematic offset from the input magnitude to mean recovered magnitude for the 
artificial stars was always much smaller than the la random errors, so we applied no systematic 
corrections to the photometry of the real objects. The errors produced by SExtractor were found 
to agree well with the errors deduced from the simulations, except for objects near the detection 
limit. In principal, there will be a surface brightness dependence since errors measured from stars 
will not be the same as those measured from less compact sources such as galaxies (e.g., Bershady 
et al. 1998). However, this effect is more important for determining the completeness corrections, 
which is not relevant to our work here, and, fTirthermore, at the faint limit the majority of our 
sources are nearly unresolved so any error in measuring the error will be negligible. 
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Since the Keck/NIRC imaging reaches fainter magnitudes than the KPNO data, they were 
processed separately, but using SExtractor in an analogous fashion. Objects were identified and 
photometric apertures were determined from the i^-band image. This information was used to 
process the other bands, with Monte Carlo simulations used to compute photometric errors. 

2.4. Optical Spectroscopy 

We observed the 1213-0017 field with the 10-m Keck II Telescope using the facility optical 
spectrograph LRIS (Oke et al. 1995) on U.T. 1999 June 15. We prepared a focal plane mask of 
l'.'5-wide slitlets to observe seven objects with TZ — K > 5. The length of each slitlet varied, with 
a minimum length of 14". We used the 150 lines/mm grating blazed at 7500 A to cover a nominal 
wavelength region from 4000-10000 A, depending on the position of the slitlets on the mask. The 
dispersion of ~4.8 A/pixel resulted in a spectral resolution of ~25 A, as measured by the FWHM 
of a few strong emission lines of serendipitous objects and the lines of the calibration arc lamp 
spectra. Sky conditions were photometric with seeing of ^.'45 at the start of the integrations. We 
obtained 3 exposures of 1800 s each with the mask oriented at a position angle of 46° east of north. 
We dithered the telescope a few arcseconds along the slit direction between each exposure. 

The slitmask data for the objects were separated into seven individual spectra and then re- 
duced using standard longslit tasks in IRAF. The data were bias-corrected using the unilluminated 
pre- and post-scan regions of the detector and flat-ficlded using internal quartz lamps obtained im- 
mediately after the observations. Sky lines were subtracted by fitting a low-order analytic function 
to each column of the images, perpendicular to the dispersion axis. We corrected each exposure 
for long-wavelength fringing by taking an average of the other 2 exposures and subtracting this 
average. The exposures were registered and combined. The most compact objects in the resulting 
mosaic are ~ 0'.'6 FWHM in the spatial direction, though they may not be point sources. A 1-d 
spectrum of each red galaxy was extracted using a 1'.'3 (6 pixel) wide aperture and was wavelength- 
calibrated using HgKr and NcAr lamps taken immediately after the observations, with an RMS of 
« 0.3 — 0.4 A in the wavelength solutions. Zeropoint shifts of < 2 A were applied to the spectra 
based on the measured wavelengths of the night sky lines. Flux calibration was performed using 
the standard star G191B2B (Massey et al. 1988; Massey & Gronwah 1990) 

3. Results 

3.1. An Excess of EROs 

Figure 1 presents the KPNO 7^-band and iC-band images of the 1213-0017 field, which cover 
an area of 11.2 arcmin"^. Figure 2 shows the field's color- magnitude diagram along with the 
no-evolution locus of local elliptical galaxies redshifted to z = 1.3, computed from multicolor 
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photometry of the Coma cluster as described in Stanford et al. (1998). If we were to account for 
passive stellar evolution, the locus would appear somewhat bluer and brighter. 

A population of objects with red {TZ — K > 5) colors are found starting at K 18 mag, 
irregularly distributed on the sky. Most of the galaxies are somewhat bluer than the Coma locus 
at z = 1.3, and the brightest of the EROs have X-band fluxes comparable to the brightest Coma 
cD galaxies. As we discuss below, there is a significant excess of red objects in this field compared 
to the surface density measured from blank sky surveys. However, before analyzing the counts of 
the reddest objects, we must first account for a few issues: (1) the effective wavelength of the TZ 
filter is significantly redder (!=a 350 A) than the standard Cousins Rc filter; (2) our data becomes 
incomplete at > 19.5, as gauged by the turnover in the number counts; and (3) the exposure time 
is not uniform over the field, with only about half the field receiving at least 50% of the maximum 
integration. (Also, there will be some added uncertainty since the abundance of the reddest objects 
has a strong dependence on the choice of color and magnitude cuts.) 

The most extensive R—K color-selected counts in a blank field of the sky come from the CADIS 
survey of Thompson et al. (1999). They find a surface density of 0.039±0.016 arcmin"^ in an area of 
154 arcmin^ for the bright EROs {R-K' > 6 and K' < 19), slightly revised to 0.047±0.009 arcmin'^ 
in their final survey area of 600 arcmin^ (D. Thompson, priv. communication). Direct comparison 
of our ERO counts with theirs is complicated as their Rcadis filter is k. 450 A bluer than our 
TZ filter. Because EROs are very red and their number counts are a strong function oi R — K 
color, the difference is significant. (The difference in IR filters, their K' versus our Ks, is much less 
important.) We computed the offset between R filters as a function of redshift using the elliptical 
galaxy SED of Coleman et al. (1980). The true SEDs of the ERO population(s) are unknown, but 
we adopt this template as a reasonable guess. At z = 0, {Rcadis — ^) is -0.08 mag; the difference 
is negative due to the fact that TZ mags are AB and Rcadis mags are Vega-based (see footnote 5). 
This difference increases nearly monotonically with redshift out to z ~ 1.5, becoming as large as 
0.4 mag. For redshifts of 0.5 — 1.7, the calculated offset is > 0.17 mag. Since the few EROs with 
known redshifts lie in this range, we add 0.17 mag to our measured TZ — K colors to compare to 
the Thompson et al. results. This is the minimum expected amount for the transformation (i.e., an 
ERO at 2; = 0.5) and therefore leads to conservative estimates for the 1213-0017 ERO overdensity. 

After accounting for this difference in filters, we find five K < 19 objects with {Rcadis— K) > 6 
in our images, four of which lie in the central 6.2 arcmin^ area which received least half of the 
maximum integration time. Since the counts over the whole mosaic are guaranteed to be incomplete, 
we consider only this central region. The counts of Thompson et al. (1999) predict 0.24 objects in a 
6.2 arcmin^ area, 16 x fewer than observed. The Poisson probability of our observing four or more 
EROs given the expected surface density is only 1.1 x 10~^, i.e., the overabundance of EROs in 
the 1213-0017 field is highly statistically significant. The probability grows to 5% if the blank sky 
counts are increased by a factor of 5.7, i.e., the 1213-0017 ERO surface density is overabundant by 
this factor at the 95% confidence level. 



-10- 



Similarly, with the less restrictive color criteria of {RcADis — K) = 5 — 6, where the blue 
cutoff is about the color of a passively evolving single-burst population at z > 1 (Figure 9), the 
Thompson ct al. (1999) data predict 2.4 objects with K < 19 while we observe 7 objects. The 
Poisson probability of finding at least this many objects is 1.2 x 10~^, and the enhancement of such 
objects in the 1213-0017 field is at least a factor of 1.4 at the 95% confidence level. 

We can combine the counts in the {RcADis — K) > 6 and {RcADis — K) = 5 — 6 bins to 
compute the enhancement of the 1213-0017 field over the blank sky counts. We compute the joint 
probability of observing the 4 and 7 galaxies in these two color bins over a range of enhancement 
factors and find the most likely excess is a factor of 4.2, with the 95% confidence range being factors 
of 2.2-7.1. Note that this calculation assumes the enhancement is the same for both color bins, 
which may not be the case. 

3.2. Morphologies of the Red Galaxies 

Figure 3 and 4 show the reduced Keck/NIRC K-hand mosaics for the fields NE and SW of 
QSO 1213-0017. Note that the images reach different limiting depths; the faintest objects in the 
NE and SW data axe Ri20.5 and 21.3 mag, respectively. We compare these data with the HST 
F^lAW images in Figure 5 to examine the observed optical and IR morphologies of the objects 
with TZ ~ K > b. The HST images are centered on the QSO so we do not have images for the 
red galaxies farther away, including unfortTinatcly Rl and R8 for which we have continuum-break 
spectroscopic redshifts (§ 3.3.2). For z = 1.3, these datasets roughly correspond to the rest-frame 
near-UV/blue (3100 - 4100 A) and far-red (8700 - 10400 A). [O II] A3727 emission lies in the 
F^lAW filter bandpass for objects at z 0.9 — 1.5. All the red objects are extended in the HST 
images meaning they cannot be galactic M dwarfs. A few of them appear to be E/SO type. The 
HST images are relatively shallow (3cr limiting surface brightness //f814VK ~ 25.7 mag arcsec"^), 
but we can still see faint emission around a number of the red galaxies. Several seem to have 
both a central spheroid and a faint component suggestive of a disk, while a few arc entirely diffuse 
emission without any compact core/bulge. None show any strong evidence for interactions. Overall, 
the detected red galaxies show remarkable consistency in their morphologies between the F^lAW 
and K-band images. We return to this issue in § 4.2.1. 

3.3. Optical Spectroscopy 

From the 2-d reduced spectra, we were able to identify spectroscopic features in five of the 
red galaxies on our slitmask, both from emission lines and continuum breaks (Table 1). To more 
accurately determine the redshifts and their errors, we used the Fourier cross-correlation technique 
of Tonry & Davis (1979) as implemented in the FXCOR task of IRAF. We linearly interpolated over 
significant residuals in the spectra resulting from imperfect cancellation of the strongest telluric 
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features. The three absorption-break galaxies were cross-correlated against the UV spectra of 
F and G dwarfs from the lUE Spectra Atlas of Wu ct al. (1991), using an average spectrum for 
each spectral type. The spectra of the emission-line galaxies were cross-correlated against the 
emission-line template of an Sc galaxy from Kinney et al. (1996). 

3.3.1. R6 and R7 — Emission-Line Galaxies 

The spectrum of R6 shows a single emission line with an observed equivalent width (EW) of 
Ri60 A (Figure 6). There is ample continuum blueward of the line so we identify it as [O II] A3727, 
placing R6 at z = 1.203, slightly in the foreground of the other redshifts in this field. Assuming 
this emission originates from H II regions of massive yoiuig stars, we infer a star formation rate of 
f« 10 /15Q Mq yr~^ using the [O II] calibration of Kennicutt (1992), comparable to the most active 
spirals in the local neighborhood (Kennicutt 1983). 

The spectrum of R7 reveals the presence of an active galactic nucleus at 2; = 1.319 (Figure 6). 
The data show strong lines of [O II] A3727 and C II] A2326 with observed equivalent widths of 
« 50 A. The C III] A1909 line is possibly detected, though this line falls near the blue end of 
the spectrum where the signal is quite low. Examination of the reduced images clearly shows a 
weak [Ne IV] A2424 emission line which lies on the red edge of the strong [O I] A6300 telluric line, 
preventing a good measurement. Strangely, no Mg II A2800 emission line is seen, even though it is 
typically at least as strong as [Ne IV] A2424 and C II] A2326 in high-z AGN (e.g.. Stern et al. 1999). 
Presumably the emission line is absorbed by gas in this galaxy though no strong Mg II A2800 
absorption feature is seen. Also, there are some possible weak absorption features in the 8800- 
9200 A range near the expected location of the higher order Balmer lines, in particular at 8896 A 
which matches the Balmer H9 transition. However, spectra become quite noisy in this region and 
the features are not overly compelling. Finally, these are no signs of an old stellar population as 
the strongest continuum breaks (B2640, B2900, D4000) and the Ca H-l-K lines are all absent. 

3.3.2. Rl, R8, and RIO — Old Continuum-Break Galaxies 

Three of the red galaxies on our slitmask — Rl, R8, and RIO — are devoid of strong emission 
lines (Figures 7), but they do show continuum breaks identifiable as the rest-frame mid-UV (2000 - 
3300 A) breaks of B2640 and B2900 at z = 1.317, 1.298, and 1.290, respectively. The breaks 
arise from metal-line blanketing on either side of the tophat-shaped feature from 2640-2720 A. 
These features, which were first suggested by Morton et al. (1977) as possibly being useful for 
redshift determinations at z > 0.75, have been used recently to identify old galaxies out to 2; ~ 1.55 
(Dunlop et al. 1996; Dey ct al. 2000). At least two of the red galaxies may also have Mg II A2800 
in absorption. These spectral features are characteristic of old stellar populations, such as Galactic 
and M31 globular clusters, M32, and the cores of local elliptical galaxies (e.g.. Rose & Deng 1999; 
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Ponder et al. 1998, and references therein). Galaxies R8 and RIO also may show the 4000 A break 
at !^ 9200 A. 

For these three galaxies, we determined the redshifts using the 4300-8000 A region of the 
spectra in the cross-correlation analysis. At longer wavelengths, the random errors increase from 
enhanced telluric OH emission, and the systematic errors also increase from imperfect removal of 
the CCD fringing pattern, due to the small number of exposures we were able to acquire. In addition 
to the choice of wavelength region, another important aspect of the redshift determination is the 
choice of stellar template. We explored template spectra ranging from spectral types FO to G8, 
using both the Fourier analysis and a direct computation of the cross-correlation between galaxy 
and template as a function of redshift. All templates produced fairly consistent redshifts using 
both methods. Figure 8 shows the rest-frame galaxy spectra with lUE stellar spectra overplotted 
for comparison. 



4. Discussion 

4.1. A Concentration of Red Galaxies at z = 1.31 

We have measured five redshifts for the red galaxies in the 1213-0017 field. Four of these lie 
near z = 1.31 with the remaining one at z = 1.203. The spectrum of QSO 1213-0017 also shows 
an Mg II absorption system at z = 1.3196 (Lanzetta et al. 1987; Stcidcl &: Sargent 1992), similar 
to the redshift of 1.319 it 0.002 for galaxy R7. Could this galaxy be responsible for the absorption? 
It is located 15" from the QSO, which is 130 h^Q kpc. Based on preliminary results from a survey 
oi z = 0.3 - 0.9 Mg II absorbers by Steidel (1993) and Steidel & Dickinson (1995), the median 
projected impact parameter between QSOs and their best candidate absorber is 45 h'^^ kpc with 
an upper bound of 80 h^Q kpc. Early results from their z > 1 survey find the absorber properties 
are not dramatically different compared to the lower redshift sample (Steidel 1995; Dickinson & 
Steidel 1996). The projected separation of R7 therefore seems too large for it to be the Mg II 
absorber, though the geometry and extent of Mg II absorption systems remain open questions (c.f. 
Steidel 1995; Charlton & Churchill 1996). In addition, our HST images show a few faint galaxies 
with projected separations from the QSO of < 80 h'^Q kpc, making them better candidates for the 
absorber. Since R7 is unlikely to be the Mg II absorber, we infer that there is at least one more 
galaxy at z = 1.3196. 

Including the Mg II absorber, there are five galaxies within an angular region of 3' and close 
together in velocity. Their redshifts have an unweighted mean of 1.309 and a standard deviation 
a = 1790 zb 570 km s~^ in the mean rest frame. ^ The full spread in redshift is 3800 km s~^ in 



'^The above results use the standard formulae. Similar results are obtained using the central location ("mean") 
and scale ("dispersion") estimators recommended by Beers ct al. (1990) for non-Gaussian velocity distributions. 
Their preferred central location measures for a sample of this size, the median and the biweight-estimated mean, give 
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the mean rest frame; for comparison, consider that 95% of the galaxies in the central 6 /i^q Mpc 
of the Coma cluster and with cz < 12,000 kms lie within ±3500 km (Postman et al. 1998). 
The considerable dispersion and range of the 1213-0017 redshifts are possibly due to velocity sub- 
structure in this field. There is a hint of this in the data: two of the galaxies (R8 and RIO) lie at 
z = 1.290-1.298 while the other three are at z = 1.317- 1.320, i.e., there is a "gap" of 2500 km s"^ 
in the mean rest frame, although there is no clear segregation of the two subsets on the sky. We 
may be seeing two separate physical entities close in redshift (e.g., filamcnts/sheets/sub-clusters) 
which may later merge into one, but this speculation lies beyond the available data. 

Interestingly, the few known z !^ 1 massive clusters show large velocity dispersions and/or signs 
of sub-structure, e.g., the 3C 324 "cluster" at z = 1.15 and 1.21 (Dickinson 1997); RX J1716.6+6708 
at z = 0.81 (Henry ct al. 1997; Gioia et al. 1999; Clowe et al. 1998); and CL 0023+0423 at z = 0.84 
(Postman et al. 1998; Lubin et al. 1998). Spectroscopically confirmed galaxies in these clusters 
have angular extents of a few arcminutes but are not strongly concentrated on the sky, which is 
also the case for the red galaxies in the 1213-0017 field. In the case of CIG J0848-F4453 at ^ = 1.27 
and the neighboring RXJ0848.9-I-4452 at z = 1.26 (Stanford et al. 1997; Rosati et al. 1999), each 
individual cluster is fairly compact but taken together these two add up to an entity with a large 
velocity range separated by only a few arcminutes. 

In addition to the galaxies with spectroscopic redshifts, we can use optical/IR colors to con- 
strain the redshifts of objects in the SW field, the only one with photometric data in all filters (see 
§ 2). The galaxies' TZJK colors are consistent with Bruzual &: Chariot (1996) models of old stellar 
populations at 2; = 1 — 2 (Figure 9), with the identification largely based on the change in TZ — J 
color. The JHK colors provide less sensitivity, but they are consistent with the same redshift 
interval (Figure 10). 

To summarize, the collection of red galaxies in the 1213-0017 field is concentrated in redshift 
at z ^ 1.31. The small number of measured redshifts prevents any definitive conclusions, but the 

redshifts do have a considerable spread. The red galaxies are spread over ~ 1.5 h'^Q Mpc on the 
sky, with a surface density about a factor of four above R — K > 5 counts and a factor of 16 above 
bright ERO counts in blank fields (§ 3.1). Although we do not know all their redshifts, it is quite 
plausible that most of the red galaxies are at ^; = 1.31, as we discuss in the next section. 

4.2. Old Stellar Populations at z ^ 1.31 

4-2.1. Origin of the Extremely Red Galaxy Colors 

With the Keck and HST findings, we now examine more closely the origin of the colors of the 
1213-0017 red galaxies and EROs. Their colors are consistent with z > 1 ellipticals (Figures 2 



z = 1.317 and 1.319, respectively. The gapper method they advocate for computing the scale gives 1840 km s 
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and 9), and we have confirmed this inference for three of them using Keck spectroscopy. For 
these galaxies there is Httle doubt the large TZ — K colors are due to old stellar populations with 
a negligible amount of dust. For the red galaxies in the SW field, their TZJHK colors are also 
consistent with unreddened z = 1 — 2 ellipticals (Figures 9 and 10); the one SW galaxy with a Keck 
spectrum, R6, supports this idea as it lies at z = 1.203, though its redshift is based on weak [O II] 
emission. (The spectrum has insufficent S/N to detect any continuum breaks, if present.) So even 
if R6 is composed of an old, red stellar population (its 7?.-band flux is only marginally brightened 
by the [O II] line), the emission line points to a small amount of simultaneous star-formation or 
weak AGN activity. In fact, R6 has quite a different appearance in optical as compared to the 
near-IR (Figure 5) which may suggest significant dust reddening in its central regions, unlike most 
of the other red galaxies (see below). 

The other galaxy with a spectroscopic redshift, R7, has an optically revealed AGN. Its 7^ — ii' 
color is far too red compared to typical Seyferts (e.g., Kotilainen & Ward 1994). Its redness might 
come from dust or from old stars, although its optical spectrum shows no signs of the latter. If the 
Balmer features in the spectrum are real, this would suggest R7 has experienced significant star 
formation in the last few ^ 10® years and would hint that the red colors are due to dust reddening. 

For the remaining 1213-0017 red galaxies, we have no direct evidence about their redshift nor 

their stellar content; all that we know is their TZ — K colors. One unlikely cause for their redness 
would be a strong emission line redshifted into the X-band; the lowest possible redshift for this 
to occur would be for Ha at z ~ 2.1 — 2.6. It seems implausible that all or even any of these 
galaxies are at this redshift given the very low surface density of strong high-redshift Ha emitters 
(e.g., Thompson et al. 1996; Beckwith et al. 1998; Teplitz et al. 1998). In fact, the more likely case 
of emission-line contamination of the broad-band magnitudes is in the 7^-band, either from Ha or 
[O II]; more than half of z = 0.5 — 1.3 field galaxies show [O II] emission (Hammer et al. 1997). In 
this case, the measured TZ — K colors would be slightly bluer than the true color of the underlying 
stellar continuum. 

We find the observed optical and near-IR morphologies are fairly regular and in good agreement 
with one another (Figure 5). Both of these tell us that dust reddening in these galaxies is spatially 
uniform on scales of < 4 h^^ kpc, as small as 0.9 h^Q kpc for the HST images. This is in 
sharp contrast to the highly wavelength-dependent appearance of the prototypical ERO HR 10: 
in the rest-frame near-UV/blue it has an inverted S-shape while in the rest-frame red it shows a 
smooth compact appearance (Dey et al. 1999). This implies HR 10 is an interacting/disturbed 
system whose central region has a large amount of dust, an interpretation augmented by its sub- 
mm emission which presumably springs from a large amount of ongoing dusty star formation. HST 
optical morphologies for the first sub-mm selected sources, which are possibly similar to HR 10, 
also show disturbed appearances (Small et al. 1998). Likewise, the rest-frame UV appearance of 
local ultraluminous IRAS galaxies (ULIRGs) are significantly different than their appearance at 
longer wavelengths (Trentham et al. 1999). 
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Therefore, we conclude the wavelength-independent appearance of the 1213-0017 red galaxies 
means they are unlikely to be dusty star-forming galaxies like HRIO or ULIRGs. Though dust 
reddening may still still play a role, most of their redness probably originates from the light of 
old stars. Dust variations could still exist on spatial scales below our resolution limits. For this 
reason, attempts to fit the broad-band SEDs are potentially misleading given the uncertainties in 
the dust geometry and distribution — such fits would be highly degenerate. Near-IR images with 
comparable spatial resolution to the HST optical data are needed to address this possibility. 

4-2.2. Rest-Frame UV Spectral Features 

Comparisons between the galaxy spectra and the lUE stellar spectra show the three continuum- 
break galaxies have real differences in their rest-frame mid-UV spectra. The spectra of RIO rises 
most rapidly to the red and has the strongest 2640 A and 2900 A breaks, while Rl's continuum 
is noticeably flatter and its breaks are the weakest, especially the 2900 A break. R8's spectrum 
is intermediate between RIO and Rl. Results from both cross correlation analysis (§ 3.3.2) and 
examination by eye find Rl is best-matched by early-type F dwarfs, R8 by mid-type F dwarfs, and 
RIO by late-type F and early-type G dwarfs (Figure 8) in this wavelength region. Furthermore, 
even when the lUE spectra are normalized only to the galaxy flux longward of 3000 A, there is still 
good agreement between the star and galaxy spectra all the way to the bluest wavelengths. This 
suggests any ongoing or recent star formation, which would be revealed by a rising blue continuum, 
is small. 

The UV continuum breaks at 2640 A and 2900 A are of special interest. This region of an old 
galaxy's SED is expected to be dominated by light from stars near the main-sequence turnoff. As 
the stellar population ages, the main-sequence turnoff mass decreases, and later-type stars produce 
the UV SED, leading to larger break amplitudes. In this fashion the rest-frame UV spectrum 
might be used to age date the stellar population, or at least the portion contributing the bulk of 
the UV flux. Spinrad et al. (1997) attempted to do this in a robust fashion for the old galaxy 
LBDS 53W091 sX z = 1.55. Given difliculties of absolute age calibration for these features and 
their uncertain metallicity dependence (Heap et al. 1998; Yi et al. 1999a), we do not attempt to 
age date the populations of the 1213-0017 galaxies. 

However, we can still study the stellar content of the spectroscopically old 1213-0017 galaxies 
by comparing the strength of these breaks with old populations in the local Universe. The spectra 
of R8 and RIO are barely amenable to this line of investigation, but the S/N of Rl's spectrum is too 
low so we exclude it. The standard indices for measuring these breaks, 2609/2660 and 2828/2921 
from Fanelli et al. (1990), are too narrow for our spectra: these use ~25 A wide bandpasses, meaning 
they include only 2-3 spectral resolution elements for our data. Instead, we use our own broader 
custom indices, W2640 and W2900. We also measure the moderate-bandwidth 2600 — 3000 index 
of Fanelli et al. (1990) to track the UV continuum color. These are all calculated in the standard 



-16- 



fashion, from the ratio of the fluxes in blue and red bandpasses and expressed in magnitudes: 

"Fa (blue) 



Index = —2.5 x log 



10 



(1) 



where Fx is the average flux density in the bandpasses. Objects with stronger breaks have larger 
values for the indices. The red bandpass of W2640 covers the entire top-hat feature redward of the 
2640 A break and ends before the Fe+Cr absorption at 2745 A; the blue bandpass has the same 
width and includes Fe II A2609 absorption and part of the broad BL A2538 feature (e.g., see Ponder 
et al. 1998). The blue bandpass of W2900 spans the Mg II A2800 and Mg I A2852 absorption, and 
the red bandpass is chosen to have the same width. Definitions and measurements of the standard 
indices and our custom ones are listed in Table 2. 

Figure 11 plots the results juxtaposed against individual main-sequence stars (Fanclli ct al. 
1992), local old spheroids (Ponder et al. 1998), and the z = 1.552 galaxy LBDS 53W091. The UV 
indices of R8 and RIO are consistent with local populations, though the W2640 strengths of R8 
and RIO differ, with the former having a weaker 2640 A break than even the local elliptical sample. 
Also, the two galaxies have diff'erent 2600 — 3000 colors, with R8 lying closer to the metal-poor 
local populations, which have bluer mid-UV colors. The similar break amplitudes of LBDS 53W091 
and RIO suggest these two galaxies are comparably old. If the analysis of 53W091's spectrum is 
correct, its inferred age is barely compatible with present cosmological parameters, even with the 
assumption of solar metallicity to reduce the inferred age (Spinrad et al. 1997). Therefore, RIO must 
have had a comparably high formation redshift, since the difference in look-back time is only 5% 
between z = 1.55 and z = 1.31, and also must have experienced little or no recent star-formation. 

For the local populations, the offset of the ellipticals from the stars and globular clusters in 
plots of the 2640 A break versus 2600 — 3000 color was noticed by Ponder ct al. (1998), who suggested 
it was due to the hot UV-cxcess (UVX) population in elliptical galaxies. The UVX is thought to 
originate from post-main sequence stars, though its specific origin remains uncertain (O'ConnclI 
1999). Interestingly, galaxy R8 does appear to be more consistent with the local ellipticals than the 
local non-UVX populations which might suggest a hot component in R8. If this could be confirmed 
for this galaxy and other high-redshift ellipticals, a detailed study of their UV spectra might uncover 
the origin of the UVX emission, since the potential UVX candidates evolve on different time scales 
(e.g.. Brown et al. 1998; Yi et al. 1999b). However, ongoing star formation, even just a small 
amount by mass, can greatly brighten a galaxy's UV flux, thereby complicating such an analysis. 

To summarize, comparison of the rest-frame UV spectra of Rl, R8, and RIO shows a real scatter 
in their stellar content even though their observed K-hand fluxes, which should roughly trace galaxy 
mass, are within a factor of two of each other. This may reflect differences in their star formation 
and/or merging histories. The UV color and break amplitudes of R8 and RIO are consistent with 
z = old stellar populations, though the data's limited S/N and spectral resolution hinders any 
more detailed comparison, e.g., studying the weaker Mg II and Mg I absorption features. At a 
redshift of 1.3, we are observing these old galaxies when the Universe was only 30% of its present 
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age. Therefore, improved data, both in S/N and spectral resolution, could probe the evolving 
composition of the oldest stellar populations from their relative youth at high redshift to their 
advanced age at z = 0. 

4-2.3. Relation to Galaxy Morphology and Implications for Formation 

The morphological results from the HST images are somewhat surprising if our suggestion 
that the 1213-0017 red galaxies contain old stars is correct. Most of the galaxies do not seem 
to be early-type ones but rather are disky or diffuse. Indeed, the one absorption-break galaxy 
with HST imaging, RIO, looks like a compact edge-on disk — this galaxy does not resemble an 
elliptical galaxy morphologically even though it does so spectroscopically. These findings argue 
that not all old EROs are necessarily ellipticals, at least in their structure. This may reflect a 
fundamental ambiguity in identifying ellipticals at z > 1 as compared to those at z = 0: are 
elliptical galaxies properly defined by their morphological or spectrophotometric properties (c.f. 
Kauffmann et al. 1996; Schade et al. 1999)? If they formed by monolithic collapse at very high 
redshift, there should be no ambiguity. But in current scenarios of hierarchical galaxy formation, 
massive early-type galaxies form late (z ~ 1 — 2 depending on the cosmology) by agglomeration 
of smaller older sub-units, so the stars of the resulting elliptical are actually much older than the 
equilibrium morphology. Simulations (Kauffmann & Chariot 1998) suggest these scenarios can still 
explain the homogeneity of the early-type cluster galaxy population out to z « 0.9 (e.g., Aragon- 
Salamanca et al. 1993; Kodama et al. 1998; Gladders et al. 1998; Stanford et al. 1998). However, 
it is at 2; > 1 when differences are predicted to become apparent, and there is little observational 
evidence in this regime due to the few number of z > 1 clusters known. Spectroscopic follow-up 
of the remaining red galaxies in the 1213-0017 field combined with their HST images should prove 
very useful in addressing this issue. 

5. Conclusions 

We have found an overdensity of EROs in the field of the z = 2.69 quasar QSO 1213-0017 
(UM 485), about a factor of 16 overdense compared to the blank field ERO surface density and at 
least of factor of 6 overdense at the 95% confidence level. The optical/IR colors of the EROs and 
numerous other red galaxies in this field are consistent with those of passively-evolving elliptical 
galaxies at z > 1. HST optical imaging shows a few of the red (JZ — K > 5) galaxies seem to have 
early- type morphologies while the remainder are either disk galaxies or diffuse objects without any 
obvious core. Their near-IR morphologies are consistent with those observed in the optical, unlike 
in the case of the prototypical ERO HR 10. This suggests that the dust extinction in these galaxies 
is either relatively smooth or that reddening by dust is not a significant cause of the colors. 

Follow-up Keck spectroscopy has measured redshifts for five red galaxies. Four lie at z « 1.31, 
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and three of these have rest-frame UV absorption-hne spectra similar to present-day elHptical 

galaxies, making this the most distant concentration of old galaxies spectroscopically confirmed to 
date. Including an Mg II absorber seen in the spectrum of the background quasar, there arc five 
spectroscopic redshifts at z ^ 1.31 with a standard deviation of 1800 km s~^ and a full range of 
3800 km s"^ in the mean rest frame. 

A number of lines of evidence possibly suggest the red galaxies in this field delineate the 
presence of a massive high-redshift cluster. The ERO surface density is enhanced above blank field 
counts, and there are five spectroscopic redshifts close together. The reddest 1213-0017 galaxies, 
three of which have the spectra of early-type galaxies, are nearly as luminous and as red as the 
brightest Coma cluster ellipticals. The red galaxies also have a large angular extent on the sky, and 
a considerable velocity spread, suggesting dynamical youthfulness. An impression of youthfulness 
is also conveyed by the roughly filamentary distribution of the red galaxies on the sky. 

Aside from the direct physical evidence, the presence of spectroscopically old EROs provides 
circumstantial support to the idea of a cluster. Elliptical galaxies in the local Universe, presumably 
the present-day counterparts of z > 1 old EROs, are highly clustered and predominantly found in 
high density regions (e.g.. Dressier 1980). Such an analogy suggests that old EROs should be found 
with others of their ilk — this idea has yet to be fully explored. Since the hallmark of rich clusters 
is the presence of old ellipticals, this z = 1.31 field may be one of the most distant rich cluster 
candidates to date. Moreover, regardless of whether this system is shown to be a genuine massive 
cluster, the concentration of EROs is likely a sign that this field is an uncommonly overdense region 
at high redshift, which warrants follow-up studies. 

Further observations are needed to develop a full physical picture of this field. A much larger 
sample of redshifts, inferred from multicolor photometric redshifts and directly measured from 
deep spectroscopy, will reveal the number of luminous galaxies at this redshift. An expanded 
spectroscopic sample will also be valuable to measure the galaxy velocity distribution and to search 
for dynamical sub-structure. X-ray observations will allow tis to look for hot intracluster gas, a sign 
of a deep gravitational potential, and also for mass sub-structure. This issue can also be addressed 
with radio observations to search for a Sunyaev-Zeldovich decrement. 

Extending galaxy cluster studies to z > 1 has special importance since current hierarchical 
formation theories predict massive galaxies assemble from smaller sub-units during this epoch (e.g., 
Kauffmann &; Chariot 1998), so we can potentially test these models by witnessing the formation 
process in situ. Measuring the high-redshift cluster abundance using wide-area surveys can provide 
strong tests of cosmological models (Eke et al. 1996; Bahcall et al. 1997). However, even before 
a large number of such clusters have been found, the few known to date are valuable sites to 
investigate the processes which drive the evolution of the oldest stellar populations, galaxies, and 
galaxy clusters. Furthermore, these clusters can also be used as testing grounds for z > 1 cluster- 
search strategies, which by necessity are derived from extrapolating the known properties of lower- 
redshift clusters. It may be that we have to abandon some common precepts formed from studying 
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local rich clusters in order to develop a complete understanding of high-redshift clusters and their 
constituent galaxies at an epoch of less than half the current age of the Universe. 
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Fig. 1.— KPNO 7^ and K-hand images of the 1213-0017 field. Each is 3^5 x 3;i. Objects with 
TZ — K > 6 are circled (Q) and those with TZ — K = 5 — 6 are marked with a square (□). The 
quasar is marked on the 7^-band image with a white cross and labelled "Q". Red galaxies with 
spectroscopic redshifts are labelled. 
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K magnitude 

Fig. 2. — Color-magnitude diagram for objects detected in the 11 arcmin^ iiT-band image shown 
in Fig. 1. The error bars for K < 19 objects are smaller than the plotting symbols. 7^-band 
magnitudes are AB-based, and K-band mags are Vega-based. There are 11 objects with Tl — K>6 
and and another 22 with TZ — K = 5 — 6. The dashed line represents a no-evolution color-magnitude 
locus for elliptical galaxies at z = 1.3 based on Coma cluster galaxies (Stanford ct al. 1998), with 
the bright end indicating the expected location of the brightest cluster galaxy. The stars (*) show 
the five objects with spectroscopic redshifts. 
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Fig. 3.— Keck K-band image of a field NE of QSO 1213-0017 (see § 2.2.2 for position) shown with 
a hnear stretch. The field is 47" x 55". The very faintest objects in the images are K 20.5 mag. 
EROs with spectroscopic redshifts are labeled. 
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Fig. 4.— Keck K-band image of a field SW of QSO 1213-0017 (see § 2.2.2 for position) shown with 
a linear stretch. The field is 48" x 48". The very faintest objects in the images are X ~ 21.3 mag. 
The one EROs in the image with a spectroscopic redshift is labeled. 




Fig. 5. — HST F814W and Keck K-hand images for most of the red galaxies in the field of 
QSO 1213-0017. The top set is from WF3, which covers the field NE of the quasar, and the 
bottom set from WF4, covering the SW field. For each set, the top row shows the greyscale HST 
image, the middle row shows contours for the HST data drawn at{^,|,|,|,|}xcj, and the bottom 
row shows the Keck image, when available. The orientation on the sky is shown at the right edge 
of each set. The images are 3'.'5 (30 h'^Q kpc) on a side and displayed with a linear stretch. The 
TZ — K color and TZ mag measured from our KPNO data are above each object. The three with 
redshifts (R6, R7, and RIO) are labeled. Most of the red galaxies have a regular appearance, and 
none are obviously interacting or disturbed. 
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Fig. 6. — Keck/LRIS spectrum of two emission-line red galaxies. Both spectra have been smoothed 
with a boxcar filter of width 7 pixels. Top: R6 shows a single strong emission line at 8210 A, which 
we assume is [O II] A3727. The feature at 8642 A is due to imperfect cancellation of a CCD fringe 
from a strong sky line. Bottom: R7 shows emission lines of C II] A2326 and [O II] A3727, and also 
a weak [Ne IV] A2424 line close to the telluric [O I] line at 5577 A. 
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Fig. 7. — Keck/LRIS spectrum of the three absorption-line red galaxies. The spectra have been 

averaged in 7 pixel bins with the formal la errors overlaid. See the text and later figures for a 
description of the redshift determination. At A > 8500 A, the noise of the spectra increases due to 
telluric OH emission lines and also imperfect removal of CCD fringes. 
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Fig. 8. — Rest-frame UV spectra of the three continuum-break galaxies, averaged in 7 pixel bins. 
The R8 spectrum is offset by -|-13 along the y-axis and the RIO spectrum by +26. Overplotted with 
the thick lines are lUE spectra of dwarfs. The lUE spectra have been normalized to the galaxy flux 
longward of 3000 A; notice the good agreement between the star and galaxy spectra all the way to 
the bluest wavelengths. This suggests the amount of any ongoing or recent star formation, which 
would be revealed by a rising blue continuum, is small. 
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Fig. 9.— UJK color-color plot of red {U - K > b) objects in the Keck/NIRC field SW of the 
quasar, the only field with IR photometry in two filters. The data point without error bars has 
errors comparable to its symbol size. Overplotted are single-burst models of Bruzual & Chariot 
(1996) passively evolving from formation rcdshifts of Zi? = 3 and 5. The open circles (o) show the 
observed model colors at z = — 2.5. The colors of the single-burst models diverge for z > 1.5. We 
also plot colors for an exponentially decaying star-forming model with an e-folding time of 6 Gyr, 
which well reproduces the colors of local spiral galaxies (e.g., Chariot & Bruzual 1991); its colors 
never become significantly red. The observed colors of the wcll-dctcctcd red galaxies suggest these 
objects lie at z = 1 — 2, corroborating the color-magnitude plot in Figure 2. The object marked 
with a star {-k) is R6, an emission-line ERO at z = 1.203. 
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Fig. 10. — JHK color-color plot for all objects in the Kcck/NIRC images SW of the quasar, the 
only field where we measured two IR colors. The filled symbols are the data with errors in both 
colors of < 0.25 mag, and the empty symbols are the lower S/N data. The triangles indicate the 
TZ — K > 5 objects from Figure 9. The median errors for each subset are plotted on one data point. 
The lines show the same theoretical models as in Figure 9. The observed colors of the well-detected 
red galaxies suggest these objects lie at z = 1 — 2. The one object in the field with a redshift, the 
emission-line ERO R6 at 2; = 1.203, is marked with a star (★) and labeled. 
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Fig. 11. — Strength of the UV-continuum breaks of red galaxies R8 and RIO compared with local 
old stellar populations and LBDS 53W091 at z = 1.55. W2640 and W2900 measure the 2640 A 
and 2900 A breaks, respectively, and 2600 — 3000 is a moderate-bandwidth UV color. Data for 
individual stars are from Fanelli et al. (1992), for Milky Way globulars from Rose &; Deng (1998), 
and for M 31 globulars and ellipticals from Ponder et al. (1998). More metal-poor globulars lie 
to the lower left and more metal-rich ones to the upper right. R8 is more similar to the more 
metal-poor local systems than RIO. The W2640 offset between the M31 globular clusters and the 
local ellipticals is believed to arise from the UV-excess population in the ellipticals. 



Table 1. Spectroscopic Redshifts in the 1213-0017 Field 



Identification 


RA (J2000) 


Dec (J2000) 


K (mag) 


TZ- K (mag)'' 


Morphology'^ 


z 


Spectral Features 


Rl 


12 15 43.2 


-00 35 39 


18.30 ±0.07 


6.42 ±0.15 




1.317 ±0.005 


B2640, Mg II, B2900 


R6 


12 15 48.1 


-00 34 37 


18.89 ±0.07 


5.96 ±0.12 


diffuse'^ 


1.203 ±0.002 


[0 II] 


R7 


12 15 50.9 


-00 34 32 


18.04 ± 0.04 


5.84 ±0.07 


compact+disk? 


1.319 ±0.002 


[0 II],CII], [NelV], C III]? 


R8 


12 15 44.8 


00 34 25 


18.94 ±0.09 


5.30 ±0.13 




1.298 ±0.004 


B2640, Mg II?, B2900, D4000? 


RIO 


12 15 52.2 


-00 34 00 


18.19 ±0.05 


6.26 ±0.08 


edge-on disk 


1.290 ±0.002 


B2640, Mg II, B2900, D4000? 


Mg II absorption^ 


12 15 49.8 


-00 34 34 








1.3196 
















1.5534 





'^T^.-band mags arc on the AB system while if -band mags are normalized to Vega (see § 2.1). A flat spectrum (/^ = constant) source has TZ—K = 1.85. 
''Morphologies from our HST F814W images. Sources with a blank line " — " were not observed. 

'^Unlike the other EROs, which have similar morphologies in the FSIAW and if -band images, R6 appears diffuse in i^814H^ but compact in the 
if -band images. See § 4.2.1 and Figure 5. 

<iProm Steidel & Sargent (1992). 
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Table 2. Rest-Prame UV Spectral Indices 



Name 


Blue Passband (A) 


Red Passband (A) 


Definition 


R8 


RIO 


2600 - 3000 

2609/2660 

2828/2921 


2470-2670 
2596-2623 
2818-2838 


2930-3130 
2647-2673 
2906-2936 


Fanelli et al. (1990) 


0.63 ± 0.09 
0.41 ± 0.34 
0.41 ±0.31 


1.22 ±0.21 
0.62 ±0.61 
0.29 ±0.36 


W2640 

W2900 


2550-2625 

2773 2873 


2647-2722 

2910 3010 


this work 


0.05 ± 0.18 

0.50 ±0.12 


0.73 ± 0.45 

0.71 ±0.23 



